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We have experimentally measured transport of superfluid, bosonic atoms in a mesoscopic system: a small 
channel connecting two large reservoirs. Starting far from equilibrium (superfluid in a single reservoir), we 
observe first resistive flow transitioning at a critical current into superflow, characterized by oscillations. We 
reproduce this full evolution with a simple electronic circuit model. We compare our fitted conductance to 
two different microscopic phenomenological models. We also show that the oscillations are consistent with 
LC oscillations as estimated by the kinetic inductance and effective capacitance in our system. Our experiment 
provides an attractive platform to begin to probe the mesoscopic transport properties of a dilute, superfluid, Bose 
gas. 


I. INTRODUCTION 

Transport phenomena in mesoscopic systems are character¬ 
ized by the importance of quantum phase coherence. In these 
systems, the length scale associated with the device is typi¬ 
cally smaller than or comparable to the inelastic mean free 
path. This can lead to a wide variety of different effects, in¬ 
cluding quantum conductance [1] and quantum persistent cur¬ 
rents in normal metal rings [2, 3]. Cold atomic gases typi¬ 
cally have mean free paths longer than the system size, and 
thus mesoscopic transport effects are crucial. For example, in 
degenerate Fermi gases, quantum conductance has been ob¬ 
served [4]. Here, we measure mesoscopic transport in a Bose- 
condensed gas and observe superflow below a critical current 
and resistive flow, possibly associated with the creation of vor¬ 
tex pairs, above this critical current. The dynamics of this 
system is well described using a simple circuit model which 
captures the essential physics. In turn, the circuit parameters 
can be used to search for a microscopic explanation or to in¬ 
form a more full modeling of system with mean-field theory. 

Our system consists of two large condensates, or reservoirs, 
connected by a channel that is long compared to the heal¬ 
ing length [5] of the condensate but small compared to the 
mean free path (Figure 1). A similar experiment with thermal 
bosons observed a ballistic (Sharvin) resistance [6]. Exper¬ 
iments in an analogous experiment [4, 7-12] with fermions 
observed the superfluid transition [8] and thermoelectric ef¬ 
fects [9]. Our system, prototypical of many mesoscopcic 
transport devices, is of theoretical interest [13, 14] because it 
may help lead to new cooling mechanisms [15, 16] and obser¬ 
vation of the superfluid fountain effect [17]. In addition, if the 
channel is in the one-dimensional regime it can be described 
with a Hamiltonian similar to a Luttinger liquid [18] and could 
violate the Wiedemann-Franz law [19]. Because of the long 
length of our channel, we expect to see different transport ef¬ 
fects compared to using either tunnel barriers or short weak 
links [20, 21]. Experiments with such junctions have observed 
quantum effects like macroscopic quantum self-trapping [22], 
the ac and dc-Josephson effects [23], and the transition from 
tunneling junctions to weak links [24]. Similar experiments 
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FIG. 1. In-situ images of a BEC of 495(16) xlO'^ ^^Na atoms in 
the dumbbell-shaped potential for trapping parameters that yield an 
equilibrium 1-D density of atoms in the channel of 790(25) pm“'. 
The atoms are initially trapped in the left reservoir. At time ? = 0 ms, 
a gate is removed, and the atoms are allowed to flow freely between 
the reservoirs. 


with weak links and tunnel junctions in rings have shown re¬ 
sistive flow in superfluids [25, 26], persistent currents [27], 
and discrete phase slips [28-30]. 

In the experiments reported here, we observe that while the 
phase difference between the condensates governs the super¬ 
fluid transport, there is also large dissipation. This dissipation 
is related to the creation of excitations in one of the reservoirs, 
an effect thought to contribute to the resistance [26] but not 
conclusively shown. Because the creation of excitations ap¬ 
pears to occur not within the channel but at the interface with 
the reservoir, the model allows us to consider our dissipation 
as a “contact resistance”. Contact resistance is a hallmark of 
mesoscopic transport: because the channel length is smaller 
than the inelastic mean free path, any dissipation must occur 
in the contacts. 


II. BRIEF EXPERIMENTAL DESCRIPTION 

We start our experiments with all the atoms contained in 
one of the two reservoirs, as shown in Fig. 1. The atoms 
are contained in this “source” reservoir by a “gate” potential. 
When the gate potential is removed, the condensate starts to 
expand through the channel. The atoms that enter the channel 
are accelerated by the change in the interaction mean field en¬ 
ergy and spray into the “drain” reservoir. They then bounce 
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FIG. 2. (a) Example 10 ms time-of-flight (TOE) images of the 

condensate after the atom number imbalance has reached equilib¬ 
rium (612 ms after opening the gate). In almost all TOE images, we 
observe vortices, primarily in the initially empty reservoir, showing 
evidence of the Feynman mechanism for vortex production [32]. (b) 
In the Feynman model, superfluid flows out of a channel into a reser¬ 
voir. When the flow rate exceeds a critical value, vortex anti-vortex 
pairs are created. 


off the walls of the reservoir and each other, causing them to 
thermalize. 

Once a superfluid EEC is established in the drain reservoir, 
we expect that the supercurrent /j between the two reservoirs 
will be related to differences between their two phases. If the 
rate is higher than the critical current of the channel, excita¬ 
tions will be created [31-43]. Such excitations remove energy 
from the flow and eventually dissipate as heat. Time-of-flight 
imaging shows vortices in the drain reservoir (Fig. 2a), and, 
while these are not the only type of excitation, these vortex 
excitations might play a key role. The existence of vortices in 
this geometry is suggestive of the process described by Feyn¬ 
man [32] and shown in Fig. 2b. There, a channel of width 
d carries atoms into a large reservoir, producing vortices at 
the corners where contact is made. Refs. [44, 45] success¬ 
fully used the Feynman model to predict the critical velocity 
of superfluid liquid helium flowing through an orifice into a 
reservoir. 

As the entire system moves toward equilibrium, the current 
drops below a critical value (the critical current) and dissipa¬ 
tion decreases dramatically. Any chemical potential imbal¬ 
ance that still exists at this time will result in number oscilla¬ 
tions between the two reservoirs. These oscillations are analo¬ 
gous to plasma oscillations in a superconducting junction [46] 
and isothermal oscillations in superfluid liquid helium trans¬ 
port experiments [47, 48]. 

Therefore, we expect to see two distinct behaviors, depend¬ 
ing on the time after release. We first expect to have a large 
amount of dissipation from the excitations that causes the 
mass imbalance to decay. After the current drops below the 
critical value, we expect the resistive flow to decrease signif¬ 
icantly and the current to oscillate. Experiments with super¬ 
fluid liquid helium in similar configurations produce qualita¬ 
tively similar behavior (for examples, see Refs. [48, 49]). 


III. EXPERIMENTAL DETAILS 

Our experimental setup consists of a EEC of 
IF = 1, m/r =-1) ^^Na atoms in an optical dipole trap 
(ODT). Our EECs are created using standard laser cooling 
techniques, followed by evaporation in first magnetic then 
optical dipole traps. The number of atoms in the condensate 
can be tuned between 10^ and 5 x 10^ atoms, corresponding 
to equilibrium chemical potentials of jUc/S ^ In x (500 Hz) 
to jUc/S xi 2n X (1000 Hz). Vertical confinement is cre¬ 
ated using a red-detuned (1064 nm) ODT in the shape 
of a sheet. This sheet has a vertical trapping frequency of 
^ 529(2) Hz [50]. The sheet also provides confinement 
in the plane of the dumbbell. This residual confinement is 
characterized by a horizontal trapping frequency of ^ 9 Hz 
along the long axis of the dumbbell trap . 

In the horizontal plane, we use a direct intensity masking 
technique [51] to create the blue-detuned (532 nm) trap in 
the shape of dumbbell. This approach uses a Gaussian beam 
passing through a photomask in the shape of our desired po¬ 
tential. The mask is imaged onto the atoms using the same 
optical system used to observe them. The full-width, half 
max of the Gaussian beam is chosen to be approximately the 
end-to-end length of the dumbbell-shaped potential. In the 
plane, the reservoirs are nearly hard-walled with a diameter 
D - 40(3) yum. They are connected by a channel with length 
I - 22(1) yum whose potential in y is complicated by imperfec¬ 
tions in the imaging process. We empirically observe that the 
apparent Thomas-Fermi width of the channel d = 6.4(2) fj.m 
is independent of both total atom number N and strength of 
the optical potential f/„,. However, the 1-D density of atoms 
in the channel, nyo, does change with both U,„ and N. 

There is also a variable-height gate potential that is used to 
block the channel in the middle. The gate potential is created 
with a blue-detuned Gaussian beam that is scanned across the 
channel at 2 kHz using an acousto-optic deflector. The com¬ 
bined height of this time-averaged gate potential and the static 
channel potential can be adjusted to be higher than the initial 
chemical potential of the source reservoir, blocking the chan¬ 
nel. 

Due to the blue-detuned dumbbell trap having local min¬ 
ima outside of the region of interest, we adiabatically transfer 
the atoms into the source reservoir from an initial red-detuned 
ODT. The atoms are held in this configuration for at least 2 s 
to equilibriate, then the gate beam is turned off suddenly and 
the system is allowed to evolve. 

At various times after opening the channel (typically rang¬ 
ing between 5 ms to < 1 s), we count the number of atoms 
in each reservoir using partial-transfer absorption imaging 
(PTAI) [52]. To count the atoms in a given reservoir accu¬ 
rately, the transfer fraction is chosen to produce images with 
maximum optical densities between 1 and 2. For t near zero, 
images that count the number of atoms in the source reser¬ 
voir will therefore have a small transfer fraction and cannot 
accurately count the small number of atoms in the drain. To 
correct for this, a second image is taken with larger transfer 
fraction to accurately count the atoms in the drain. Therefore, 
for each time observed, a pair of images is taken to determine 
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the atom number in each reservoir. In general, the resulting 
atom numbers from three or four parrs of images are averaged 
to determine the atom number imbalance at a given discharge 
time. 

Because determining the 2D density from the optical den¬ 
sity requires dividing by the transfer fraction, atomic densities 
determined from images with smaller transfer fractions will 
have larger noise. This effect can be seen in Fig. 1 . The im¬ 
ages used to determine the atomic density in the source reser¬ 
voir have a smaller transfer fraction and larger relative noise 
than the images used to determine the atomic density in the 
drain. 

For these experiments, the temperature of our condensate 
is ^ 100 nK. This temperature is determined by time-of-flight 
absorption imaging in a direction parallel to the plane of the 
dumbbell. In this direction, the optical density of the thermal 
component is sufficiently large for detection. In the plane of 
the dumbbell, the thermal component is too sparse to be ef¬ 
fectively imaged. We estimate the critical temperature of our 
condensate to be ^ 500 nK, and thus > 95 % of the atoms 
are in the condensed state. Therefore, we expect little contri¬ 
bution to the bulk transport from the remaining 5 % thermal 
atoms. Moreover, the thermal cloud exists in regions outside 
the dumbbell (but still confined by the sheet potential), so we 
expect reasonable thermal contact between the two reservoirs. 

The process of discharge represents a conversion of energy 
from chemical potential to both kinetic energy (in the form of 
collective excitations) and eventually to thermal energy. We 
can estimate the maximum temperature increase of the system 
by considering the total energy of condensate E ~ Nfj.. The 
difference in the total energies of the initial state (fi/ks ~ 30 
to 80 nK, depending on parameters) and final state (/r/kg « 
20 to 50 nK) per particle represents the maximum increase in 
temperature. For the parameters of our system, we expect this 
increase to be 10 to 20 nK. 


IV. DATA AND CIRCUIT MODEL 

Figure 3 shows the atom number imbalance between the 
two reservoirs as a function of time after the gate potential is 
switched off. Here we define Ne, the equilibrium number of 
atoms in either reservoir. We further define AN as the number 
imbalance, making Ne + AN (Ne - AN) the number of atoms in 
the source (drain) reservoir. The plotted value, n = ANjlNe, 
is the normalized atom number imbalance, and can vary be¬ 
tween -1 and 1 [where n > 0 (n < 0) represents more atoms 
in the source (drain) reservoir]. The evolution for three differ¬ 
ent niD with 472(22) x 10^ atoms are shown. As predicted, 
the atom number imbalance undergoes a short-time decay fol¬ 
lowed by an oscillation. (No statistically significant conclu¬ 
sion can be drawn with respect to the timescale of the decay.) 
We note that nio (not shown) reaches an equilibrium quickly, 
typically < 20 ms. 

To fit this behavior, we model our system as a circuit that 
captures the essential physics described above. Specifically, 
we consider the circuit in the inset of Fig. 3, which is a ca¬ 
pacitor C [6] that discharges through an inductor L connected 



FIG. 3. Normalized atom number imbalance between the two reser¬ 
voirs vs. time, for a total atom number of 472(22) x lOk The three 
plots shown are for different one-dimensional densities of atoms 
in the channel uid: red circles, 790(25) pm^'; green triangles, 
665(16) yum“‘; and blue, inverted triangles 599(17) pm“*. For clar¬ 
ity, decay measurements are artificially offset vertically by 0.5. The 
solid curves are fits to the expected dynamics from the circuit shown 
in the inset (see text). 


in series to a resistance-shunted weak link (Josephson junc¬ 
tion). The capacitor represents energy stored in the chemical 
potential difference Afi between the two reservoirs, while the 
inductor represents kinetic energy stored in the flow of atoms, 
both outside and inside the channel. (C is the only parameter 
we calculate a priori, see Sec. IV A.) The weak link sets the 
critical current of the superfluid, and the resistor in parallel 
allows additional current to flow, but with dissipation. 

There are three dynamical variables in this circuit; the num¬ 
ber imbalance on the capacitor AN (defined above), the super¬ 
fluid phase difference across the weak link y, and the number 
current I. The corresponding differential equations are 


d(AN) 

dt 


( 1 ) 


n^^v^R[i-iefir)] 


\AN 


\AN J 


— “ 

— +R(I-Ief(y)) 


( 2 ) 

( 3 ) 


where /j = hfiy) is the current-phase relationship of the 
weak link, V is the voltage across the resistor and weak link, 
and h - hjln. Equation 1 defines the current, Eq. 2 is the 
AC Josephson law, and Eq. 3 is Kirchhoff’s law for voltage 
around the full circuit. These equations are integrated numer¬ 
ically with three independent parameters t = RC, uP' - 11LC, 
and Ic, which are determined by fitting to the data. The 
curves shown in Eig. 3 are the best fits for these data assum¬ 
ing f(y) - siny. (As long as /(y) is 2;7r-periodic, differ¬ 
ent current-phase relationships do not significantly alter the 
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FIG. 4. Calculated reservoir chemical potential fi (green points) 
vs. atom number N. The magneta line shows the best fit power law. 
Inset: measured power required in the gate beam to trap all the atoms 
in the source with Thomas-Fermi calculated chemical potential ji. 
The red, dashed line shows a linear fit and confirms the expected 
scaling. 

curves.) By calculating C, we are able to extract the values of 
R and L from the fitted parameters in our model, t and u. 

A. The capacitance 

We calculate the chemical capacitance following the meth¬ 
ods used in Ref. [6] using our knowledge of the reservoir 
potential (see Appendix A 1). We calculate numerically, us¬ 
ing the Thomas-Fermi approximation, how the chemical po¬ 
tential changes as a function of atom number for a fixed 
Umiti = 2nxl .8 kHz. The result is shown by the green points 
in Fig. 4. The calculated is best fit by a power law of the 
form ^ - a + pN'^, where y - 0.52(1), ajh - 45(1) Hz, and 
pjh - 1.60(2) Hz. While y is roughly independent of t/„, the 
constants a and p depend strongly on Um- 

By measuring the chemical potential relative to the height 
of our gate beam, we can experimentally verify this scaling 
behavior. In particular, we measured the power of the gate 
beam that is required to prevent any atoms from spilling from 
the source reservoir into the drain for t/„, = 1.8(4) kHz. The 
data is shown in the inset of Fig. 4, plotted versus the calcu¬ 
lated Thomas-Fermi chemical potential. The scaling should 
be linear, with an offset in fi that is roughly given by the offset 
in the channel potential (as the gate is placed across the chan¬ 
nel and thus adds to it). The fitted offset value of 300(100) Hz 
agrees with our knowledge of the channel potential (see Ap¬ 
pendix A 2). 

The chemical capacitance of the system can now be cal¬ 
culated from the difference in chemical potential between the 
two reservoirs: 

Afi^p {{Ne + ANY - (Ne - ANY) (4) 

^ {PNJ) 2r^ = 2r^^ AA (5) 


With y — 0.52, the linear approximation made here repre¬ 
sents less than a 10% error over a number imbalance up to 
approximately ANfNe - +0.80. With this approximation, the 
chemical capacitance of our system is then given by 

. ( 6 ) 

Afi 2y(iUe - a) 

B. The weak link 

The phase between the two condensates is approximated as 
being well defined in our model, as the variable y is the phase 
difference across the weak link. The 27r periodicity of f{y) 
causes the supercurrent to oscillate when the chemical poten¬ 
tial difference across the weak link V is large. This is the ac- 
Josephson effect. Thus, the decay portion of the dynamics is 
similar to the self-trapped regime observed in Refs. [22, 23]. 
However, due to the inductance, the high frequency (initially, 
about 500 Hz) oscillating current does not appear on the ca¬ 
pacitor, and the observed total current does not contain visible 
Josephson oscillations, even in the model. 

Given the nature of the discharge, one might wonder if ap¬ 
proximating the superfluid flow as being through a weak link 
is accurate. For example, at times before a superfluid is es¬ 
tablished in the drain reservoir, the phase difference should be 
undefined and the flow should be completely resistive. The 
model works at these early times because the modeled super¬ 
fluid flow is oscillating rapidly and averages to zero. At later 
times when the oscillation is dominant, the phase profile in the 
drain is complicated by excitations, yet the model effectively 
averages the phase over the entire reservoir. This two-mode 
approximation [53] works because the local phase fluctuations 
in the drain due to excitations fluctuate at timescales that are 
smaller than the oscillation frequency. A reasonable estimate 
of this timescale is the time for a phonon to traverse the reser¬ 
voir, i.e., the diameter of the reservoir divided by the speed of 
sound. For our system, this ranges from 10 to 20 ms, depend¬ 
ing on the parameters. Thus, over the period of oscillation, the 
fluctuations of the phase in the drain average out. 


C. The conductance 

Fig. 5 shows our extracted conductances G - \ jR. In meso¬ 
scopic transport, the conductance is generally determined by 
the number of available single-particle transport modes [54]. 
A fermionic system with a small number of transport modes 
(i.e., transverse single-particle modes that have an energy cut¬ 
off less than the chemical potential in either reservoir) displays 
quantum conductance, with each channel contributing a con¬ 
ductance of /i ', where h is Planck’s constant. Our particles 
are condensed bosons, and the absence of the Pauli exclusion 
principle implies no restriction on the conductance of an indi¬ 
vidual channel. 

We calculate the number of single particle transport modes 
that are available in the channel. As described in Ap¬ 
pendix A 2, the exact form the channel potential is unknown, 
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FIG. 5. Measured conductance G vs. (a) the one-dimensional den¬ 
sity of atoms in the channel, wid, and (b) the height of the optical po¬ 
tential, The solid black line is a linear fit, and the colored curves 
are the fit to the Feynman conductance (see text). The experimental 
data are grouped by total atom number (violet squares 472(22) x 10^; 
cyan, inverted triangle 331(11) x 10^; green triangles 229(9) x 10^; 
red circles 125(6) x 10^). The violet squares correspond to the data 
shown in Fig. 3. 


but is best described with V oc 3 ;“^ with a non-negligible offset 
that has a component that is proportional to Um and a compo¬ 
nent that is constant. Because the exact details of the channel 
remain unknown, we instead make a conservative estimate of 
the number of open transport channels by neglecting the off¬ 
sets and considering the channel as being harmonic in the z 
direction and a square well with width d in the y direction. In 
this case, the energies of the single particle states with zero 
transverse momentum are given by 

n^tp-n^. I 1 \ 

Eiriy, n^) = ''' 2 I 

where Uy and are the quantum numbers of the 2-D harmonic 
oscillator and m is the mass of an atom. The number of trans¬ 
port modes is then given by the number of combinations of n,. 
and that satisfy 


E(ny,nP)<fi. ( 8 ) 

Our system thus has between 3 and 11 transport channels, de¬ 
pending on and Um, yet we observe conductance of up to 

2000 h-P 

As shown in Fig. 5a, we observe that G is directly propor¬ 
tional to Hid. Because I - GV =» GA// during the decay, this 
dependence implies that the average velocity of the atoms in 
the channel is directly proportional to AyU. Alternatively, be¬ 
cause Gfh is unitless, Gjh - anjo implies the existence of a 
constant length scale a = 3.7(4) jum in the system. 

Ref. [26] suggested a simple model of phase-slip dominated 
conductance. In this model, the conductance is attributed 
to the creation of excitations that carry nizjf atoms, where 
^ is the condensate healing length and thus is 

the relevant length scale for excitations. The resulting con¬ 
ductance is Gpsjh - 2«io^. We find a good fit with a single 
scaling parameter, G = apsGps where aps = 3.9(4), but only 


when using the local ^ in the reservoirs (as opposed to the 
local ^ in the channel). 


1. The Feynman model 


In the model that Feynman describes, superfluid flows at a 
constant rate through a channel into an infinite reservoir, as 
shown in Fig. 2b. Above some critical velocity, the fluid can 
no longer sustain superflow, and vortices will be produced in 
pairs at either side of the channel, dissipating energy. Our 
experimental setup provides a unique opportunity to observe 
a system similar to that envisioned by Feynman. 

This model calculates the rate of vortex production jp as 

v^m I^m 
2nn 2 ;rfinjp 


vortex pairs per second, where v is the flow velocity. In the 
reservoir, we can also estimate the energy of a pair of vortices 
separated by a distance w as [32] 


Ep- 


7m2DlP 


In 



( 10 ) 


where n 2 D is the 2-D superfluid density. This calculation is 
for a vortex anti-vortex pair in a homogeneous BEC. To ac¬ 
count for our finite geometry, we can use a method of images 
similar to that used in electro-magnetism [55] to numerically 
calculate the energy of a pair of vortices close to the wall of 
a circularly contained BEC. This calculation involves intro¬ 
ducing image vortices outside of the reservoir such that the 
velocity field satisfies the boundary condition that it is tangent 
to the edges of the reservoir. Calculating the energy of the 
BEC with the resulting velocity field gives a correction factor 
of A- ai 1.7 to the energy in Eq. 10. Because the vortex pairs 
will be created at either side of the channel, we setw-d in 
this calculation, where d is the width of the channel. 

Erom Eqs. 9 and 10, the rate of energy dissipation in the 
system will be 


= = ( 11 ) 

Equating this power with the power dissipated by a current 
through a resistor, P - I^Rp, we can now define the Eeynman 
resistance in our system as 


Rp = 


1 

Gp 


hnxD 

^ 

^”lD 


In 



( 12 ) 


The critical current can be calculated by equating the power 
in the flow of atoms in the channel, 

to the power dissipated by vortices, Eq. 11 [32]. Eor sim¬ 
plicity, we use the equilibrium value of the chemical potential 
fie to determine Eor fitting to the experimental data, we use 
the experimentally determined densities njo and n 2 D- Because 
Mid ~ dn 2 D, Gp is approximately proportional to iiiu- 
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FIG. 6. A comparison between experimentally observed and the¬ 
oretical values for both (a) the period of LC oscillation and (b) the 
critical velocity for vortex production in the system. For both plots, 
the data points are the values for each decay measurement (run), the 
solid line is a weighted mean of these values with a standard devi¬ 
ation represented by the shaded region. The experimental data are 
grouped by total atom number and delineated by color and symbol 
as in Fig. 5. 


We fit this theory to the data through G - apGp, where ap 
is a fit parameter. The results of this fit are shown in Fig. 5b. 
The model somewhat captures the trend, and the best ht pa¬ 
rameter is ap - 0.47(10). The discrepancy from ap - \ may 
be due to the approximate nature of the Feynman model, the 
nature of our channel and/or the fact that there are other types 
of excitations not considered. 


D. The inductance and critical current 

After the current drops below the critical current, the dis¬ 
sipation drops signihcantly, and we expect oscillations, de¬ 
scribed as plasma oscillations in Ref. [16]. Such oscillations 
represent energy coherently oscillating between kinetic en¬ 
ergy (atoms moving through the channel) and chemical po¬ 
tential differences. (Similar oscillations in liquid helium ex¬ 
periments are known as isothermal oscillations [47] or plasma 
oscillations in Josephson junctions.) Equivalently, it can be 
considered the hrst normal mode of our trap. Such oscilla¬ 
tions occur with statistical signihcance in roughly half of our 
decay measurements. All oscillations shown in Fig. 3 are sig- 
nihcant. 

A good estimate of the kinetic inductance in our sys¬ 
tem is given by the kinetic inductance of the channel Lc - 
mllnip, [6]. Using this and our calculation of the capacitance, 
we can estimate the expected oscillation frequency for our os¬ 
cillator. Figure 6a shows a comparison of the experimentally 
measured frequency to that predicted by 1 / ^JLcC. We hnd the 
experimentally measured frequencies are, on average, ^ 15% 
lower than predicted. 

Figure 6b shows a comparison of the measured critical cur¬ 
rent, extracted from the amplitude of the oscillation, and the 
critical current predicted by the Feynman model. We consider 
only decay measurements where the oscillation is statistically 
signihcant. (The amplitude of oscillation is a complicated 
function of the critical current, resistance, and capacitance. 
Therefore, lack of a clear oscillation does not indicate a zero 


critical current.) Our average measured critical current is ap¬ 
proximately a factor of 5 below the predicted value. This sug¬ 
gests that other excitations with lower critical velocities may 
be playing a role. 


V. CONCLUSION 

Our experiment provides an attractive platform to probe the 
mesoscopic transport properties of a dilute, superfluid, Bose 
gas through a small channel. The mechanism of vortex pro¬ 
duction as described by Feynman [32], scaled by a factor of 
2, predicts the general trend in our data (Fig. 5). Combined 
with the direct observation of vortices (Fig. 2), our experi¬ 
mental data suggests that the Feynman mechanism for vortex 
production plays a role in determining the conductance. To 
conclusively show the relevance of the Feynman mechanism 
in similar mesoscopic cold-atom experiments, future experi¬ 
ments should use an initial condition where both reservoirs are 
partially filled and current bias the system by contracting one 
reservoir and expanding the second. Lastly, when the current 
drops below the critical value, we observe plasma oscillations 
that are not visibly damped. 

Our experiment sets the stage for a number of other exper¬ 
iments. First, with tighter channel potentials or fewer atoms, 
the number of open transport channels could be further re¬ 
duced. By controlling the initial imbalance, it might then be¬ 
come possible to see changes in the conductance as transport 
channels are opened or closed. Moreover, this geometry could 
be used in studies of unique cooling mechanisms [15, 16], the 
superfluid fountain effect [17], or possibly seeing quantized 
conductance with bosons [56]. Because this experimental ar¬ 
rangement is capable of producing hundreds of vortex pairs, it 
could also prove useful in studying quantum turbulence [57] 
and perhaps the emergence of states like the Onsager vor¬ 
tex [58]. 
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Appendix A: Modeling the dumbbell potential 

Here, we model the potentials of the reservoir and channel 
to determine simple functional forms and best-fit parameters. 
We note, however, that in the present paper we do not rely 
on these best fit parameters in our understanding of the con¬ 
ductance or the Eeynman model, as we have explicitly written 
our models to use the measured riip, of the channel, n 2 D and 
the apparent Thomas-Eermi width d of the channel. We in¬ 
clude these models to better understand subtle effects (such as 
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the effect of the residual transverse confinement on the capac¬ 
itance of the system) and for possible future modeling of our 
experiment. 


ro = 27(2) jim. Given the numerical aperture of the imaging 
stack, the expected values of w is « 3 fim. The anomalously 
large value is most likely due to imaging aberrations, as de¬ 
scribed in detail in the next the subsection. 


1. The reservoirs 


2. The channel 


As described in Sec. Ill, the reservoir potential is made of 
three components: the harmonic confinement in the z direc¬ 
tion, the residual transverse harmonic confinement in the jc-y 
plane, and the near square well potential in the x-y plane cre¬ 
ated by the dumbbell potential itself The residual transverse 
harmonic confinement has its center in between the two reser¬ 
voirs of the dumbbell. As can be gleaned from Fig. 1, this 
residual harmonic confinement is a perturbation on the over¬ 
all potential, as it does not drastically alter the measured 2D 
density of atoms along the long axis of the channel. 

As for the box portion of the potential, diffraction and 
imaging imperfections cause the hard-walled nature to be 
smoothed out. If we were to approximate the aberrated point 
spread function of our imaging system as a Gaussian with 1 /e^ 
radius w, then the resulting box portion of the reservoir po¬ 
tential will be given by the convolution of that point spread 
function with the optical mask used to generate the potential. 
The resulting form of the reservoir potential is then 


2 2 ^ 


1 2 2 
-moj.z H-—- 
2 ^ 2 


1 + erf 




(Al) 


where r - yjix - ycf' + {x - y^)^ is the radial coordinate rela¬ 
tive to the center of the reservoir, Xc (jc) is the x (y) coordinate 
of the center of the reservoir, ro is the radius of the reservoir, 
erf is the error function, and w, is the trapping frequency in 
the /th direction. 

We use this form of the potential along with the known 
number of atoms to calculate an expected 2D density. We then 
fit, using all available equilibrium densities of both reservoirs, 
for the parameters cOy, ro, and w. The best fit values are 
toJ27T = 9.1(9) Hz, tOy/27T = 9.4(6) Hz, w = 12(2) pm, and 


If the imaging process were perfect with infinite resolu¬ 
tion (i.e, no diffraction), we expect that the channel potential 
would be given by a square well with a width d ^ \A pm. 
However, even in the absence of aberrations, our imaging sys¬ 
tem would produce an approximate square well with walls that 
changed from zero to the maximum height t/,„ over a length 
scale ss 3 pm. The potential is further complicated by the pres¬ 
ence of optical aberrations, including both spherical aberra¬ 
tion and astigmatism. As a result, a region of sightly depleted 
density appears in the channel along the long axis nearly in 
the center. This “ridge” is visible if one looks carefully at 
Fig. 1. These aberrations make effective modeling of the po¬ 
tential a priori virtually impossible, as it is unclear whether 
these aberrations are present in both imaging of the potential 
and subsequent imaging of the atoms, and, if so, in what rela¬ 
tive quantities. 

Instead, we choose to model the channel potential phe¬ 
nomenologically, using observables such as the apparent 
Thomas-Fermi width, integrated ID density nm, and cross 
sectional profile (density vs. y). We find the potential is best 
described by V oc y^ potential, with the bottom of the potential 
having an offset given by bUm + Vq, where b - 0.15(2) and 
Valh - 223(30) Hz. Here, b represents the contribution to 
the offset due to imaging aberrations and Vo represents a con¬ 
stant background potential, most likely due to a localized high 
point in the potential generated by the sheet beam. Note that 
to accurately reproduce the data, one must take into account 
the 2D-3D crossover: if p - Viy = 0) < bo)^, we use only the 
ground state of the harmonic oscillator in the z direction; if 
p - Viy = 0) > bo)^, we use the Thomas-Fermi solution in the 
z direction. 
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